###This R script replicates the models from Table A10 from the Online Appendix
### Load packages
require(foreign)
require(sandwich)
require(lmtest)
require(MASS)
require(Hmisc)
require(stargazer)
require(sciplot)
require(lme4)

### Clear space
rm(list = ls())


# define function to combine multe model estimates (this incorporates uncertatiny from the latent variables, see the supplementary appendix section F for more details about this function)
milm <- function(fml, midata){
  xx <- terms(as.formula(fml))
  lms <- matrix(data=NA, nrow=(length(attr(xx, "term.labels")) + 1), ncol=length(midata))
  ses <- matrix(data=NA, nrow=(length(attr(xx, "term.labels")) + 1), ncol=length(midata))
  vcovs <- list()
  for(i in 1:length(midata)){
    tmp <- lm(formula=as.formula(fml), data=midata[[i]])
    lms[,i] <- tmp$coefficients
    ses[,i] <- sqrt(diag(vcov(tmp)))
    vcovs[[i]] <- vcov(tmp)
  }
  par.est <- apply(lms, 1, mean)
  se.within <- apply(ses, 1, mean)
  se.between <- apply(lms, 1, var)
  se.est <- sqrt(se.within^2 + se.between*(1 + (1/length(midata))))
  list("terms"=names(tmp$coefficients), "beta" = par.est, "SE"=se.est, "vcovs"=vcovs,"coefs" = lms)    
}

### Load data. 
setwd("C:/Users/Administrator/Dropbox/Papers/Spatial repression and democracy/Data/Data Analysis")

dat <- read.dta("MainData.dta")

attach(dat)
datfh <- na.omit(data.frame(theta_mean, theta_sd, l_fariss, l_fariss_sd, SLreg_fariss_t1, SLreg_fariss_sd, v2cltort_sd, v2clkill_sd,
                            l_polity, l_polyarchy, l_polyarchy_sd, SLreg_polyarchy_t1, SLreg_polyarchy_sd, l_uds, l_boix, l_cgv, l_gwf_democracy, 
                            SLreg_polity, SLreg_gwf, SLreg_polyarchy_t1, SLreg_boix_t1, SLreg_cgv_t1,
                            SLreg_clkill_sd, SLreg_clkill_t1, SLreg_cltort_sd, SLreg_cltort_t1, l_v2clkill_sd, l_v2clkill, l_cltort_sd, l_cltort, v2cltort, v2clkill,
                            l_uds2, l_uds_sd, SLreg_uds_t1, SLreg_uds_sd, l_gdp_gro, l_lngdppc, l_lnpop, civil_war, international_war, l_resdep,
                            ccode_fe1, ccode_fe2, ccode_fe3, ccode_fe4, ccode_fe5, ccode_fe6, ccode_fe7, ccode_fe8, ccode_fe9, 
                            ccode_fe10, ccode_fe11, ccode_fe12, ccode_fe13, ccode_fe14, ccode_fe15, ccode_fe16, ccode_fe17, 
                            ccode_fe18, ccode_fe19, ccode_fe20, ccode_fe21, ccode_fe22, ccode_fe23, ccode_fe24, ccode_fe25, 
                            ccode_fe26, ccode_fe27, ccode_fe28, ccode_fe29, ccode_fe30, ccode_fe31, ccode_fe32, ccode_fe33, 
                            ccode_fe34, ccode_fe35, ccode_fe36, ccode_fe37, ccode_fe38, ccode_fe39, ccode_fe40, ccode_fe41, 
                            ccode_fe42, ccode_fe43, ccode_fe44, ccode_fe45, ccode_fe46, ccode_fe47, ccode_fe48, ccode_fe49, 
                            ccode_fe50, ccode_fe51, ccode_fe52, ccode_fe53, ccode_fe54, ccode_fe55, ccode_fe56, ccode_fe57,
                            ccode_fe58, ccode_fe59, ccode_fe60, ccode_fe61, ccode_fe62, ccode_fe63, ccode_fe64, ccode_fe65, 
                            ccode_fe66, ccode_fe67, ccode_fe68, ccode_fe69, ccode_fe70, ccode_fe71, ccode_fe72, ccode_fe73, 
                            ccode_fe74, ccode_fe75, ccode_fe76, ccode_fe77, ccode_fe78, ccode_fe79, ccode_fe80, ccode_fe81, 
                            ccode_fe82, ccode_fe83, ccode_fe84, ccode_fe85, ccode_fe86, ccode_fe87, ccode_fe88, ccode_fe89, 
                            ccode_fe90, ccode_fe91, ccode_fe92, ccode_fe93, ccode_fe94, ccode_fe95, ccode_fe96, ccode_fe97, 
                            ccode_fe98, ccode_fe99, ccode_fe100, ccode_fe101, ccode_fe102, ccode_fe103, ccode_fe104, ccode_fe105, 
                            ccode_fe106, ccode_fe107, ccode_fe108, ccode_fe109, ccode_fe110, ccode_fe111, ccode_fe112, ccode_fe113, 
                            ccode_fe114, ccode_fe115, ccode_fe116, ccode_fe117, ccode_fe118, ccode_fe119, ccode_fe120, ccode_fe121, 
                            ccode_fe122, ccode_fe123, ccode_fe124, ccode_fe125, ccode_fe126, ccode_fe127, ccode_fe128, ccode_fe129, 
                            ccode_fe130, ccode_fe131, ccode_fe132, ccode_fe133, ccode_fe134, ccode_fe135, ccode_fe136, ccode_fe137, 
                            ccode_fe138, y_fe1, y_fe2, y_fe3, y_fe4, y_fe5, y_fe6, 
                            y_fe7, y_fe8, y_fe9, y_fe10, y_fe11, y_fe12, y_fe13, y_fe14, y_fe15, y_fe16, y_fe17, y_fe18, y_fe19, y_fe20, 
                            y_fe21, y_fe22, y_fe23, y_fe24, y_fe25, y_fe26, y_fe27, y_fe28, y_fe29, y_fe30, y_fe31, y_fe32, y_fe33, y_fe34, 
                            y_fe35, y_fe36, y_fe37, y_fe38, y_fe39, y_fe40, y_fe41, y_fe42, y_fe43, y_fe44, y_fe45, y_fe46, y_fe47, y_fe48, 
                            y_fe49, y_fe50, y_fe51, y_fe52, y_fe53, y_fe54, y_fe55, y_fe56, y_fe57, y_fe58, y_fe59, y_fe60, y_fe61))

dim(datfh)


### Setup data to run models where we account for uncertainty
set.seed(32261991)
newdata <- list()
for(i in 1:1000){
  newdata[[i]] <- datfh
  newdata[[i]]$draw <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$theta_mean, sd=datfh$theta_sd)
  newdata[[i]]$draw.lag <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_fariss, sd=datfh$l_fariss_sd)
  newdata[[i]]$draw.lag.SLfariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_fariss_t1, sd=datfh$SLreg_fariss_sd)
  newdata[[i]]$draw.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2clkill, sd=datfh$v2clkill_sd)
  newdata[[i]]$draw.lag.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_v2clkill, sd=datfh$l_v2clkill_sd)
  newdata[[i]]$draw.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2cltort, sd=datfh$v2cltort_sd)
  newdata[[i]]$draw.lag.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_cltort, sd=datfh$l_cltort_sd)
  newdata[[i]]$draw.lag.SLkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_clkill_t1, sd=datfh$SLreg_clkill_sd)
  newdata[[i]]$draw.lag.SLtort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_cltort_t1, sd=datfh$SLreg_cltort_sd)
  newdata[[i]]$draw.lag.SLpolyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_polyarchy_t1, sd=datfh$SLreg_polyarchy_sd)
  newdata[[i]]$draw.lag.polyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_polyarchy, sd=datfh$l_polyarchy_sd)
}


##Model 1 - Torture (NS)
FML_NS <-"draw.cltort ~ draw.lag.cltort + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep"

model_ns <- milm(fml=FML_NS, midata=newdata)
df <- 5870

# create object for printing during the loop and then print the results to screen
model1<- data.frame(model_ns$terms, model_ns$beta, model_ns$SE, model_ns$beta/model_ns$SE)
colnames(model1) <- c("Variable","Est. Coef.","SE","z-score")

print(paste("tortings"))
print(model1)

#Model 3 - Torture Spatial Autoregressive (SAR)
FML_SAR <-"draw.cltort ~ draw.lag.cltort + draw.lag.SLtort + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"

model_sar <- milm(fml=FML_SAR, midata=newdata)

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("tortings SAR"))
print(model3)


#Model 4 - Torture Spatial Lag of X (SLX)

FML_SLX <-"draw.cltort ~ draw.lag.cltort + draw.lag.SLpolyarchy + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_slx <- milm(fml=FML_SLX, midata=newdata)

model4<- data.frame(model_slx$terms, model_slx$beta, model_slx$SE, model_slx$beta/model_slx$SE)
colnames(model4) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Kllings SLX"))
print(model4)


#Model 5 - Torture Spatial Durbin Model (SDM)
FML_SDM <-"draw.cltort ~ draw.lag.cltort + draw.lag.SLtort + draw.lag.SLpolyarchy + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_sdm <- milm(fml=FML_SDM, midata=newdata)

model5<- data.frame(model_sdm$terms, model_sdm$beta, model_sdm$SE, model_sdm$beta/model_sdm$SE)
colnames(model5) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("tortings SDM"))
print(model5)

### Tables A12 and A13 ###


##Model 1 - Killings NS
FML_NS <-"draw.clkill ~ draw.lag.clkill + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep"

model_ns <- milm(fml=FML_NS, midata=newdata)
df <- 5870

# create object for printing during the loop and then print the results to screen
model1<- data.frame(model_ns$terms, model_ns$beta, model_ns$SE, model_ns$beta/model_ns$SE)
colnames(model1) <- c("Variable","Est. Coef.","SE","t-statistic")

print(paste("Killings"))
print(model1)



#Model 3 - Killings Spatial Autoregressive (SAR)
FML_SAR <-"draw.clkill ~ draw.lag.clkill + draw.lag.SLkill + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"

model_sar <- milm(fml=FML_SAR, midata=newdata)

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Killings SAR"))
print(model3)


#Model 4 - Killings Spatial Lag of X (SLX)

FML_SLX <-"draw.clkill ~ draw.lag.clkill + draw.lag.SLpolyarchy + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_slx <- milm(fml=FML_SLX, midata=newdata)

model4<- data.frame(model_slx$terms, model_slx$beta, model_slx$SE, model_slx$beta/model_slx$SE)
colnames(model4) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Kllings SLX"))
print(model4)


#Model 5 - Killings Spatial Durbin Model (SDM)
FML_SDM <-"draw.clkill ~ draw.lag.clkill + draw.lag.SLkill + draw.lag.SLpolyarchy + draw.lag.polyarchy +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_sdm <- milm(fml=FML_SDM, midata=newdata)

model5<- data.frame(model_sdm$terms, model_sdm$beta, model_sdm$SE, model_sdm$beta/model_sdm$SE)
colnames(model5) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Killings SDM"))
print(model5)

### Section A4 - Regional conflict ###

### Tables A17 and A18 ###

##Model 1 - Non-spatial (NS)
FML_NS <-"draw ~ draw.lag + draw.lag.polyarchy + SLreg_civilwar+
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep"

model_ns <- milm(fml=FML_NS, midata=newdata)
df <- 5870

# create object for printing during the loop and then print the results to screen
model1<- data.frame(model_ns$terms, model_ns$beta, model_ns$SE, model_ns$beta/model_ns$SE)
colnames(model1) <- c("Variable","Est. Coef.","SE","t-statistic")

print(paste("Model Output"))
print(model1)


#Model 3 - Regional Civil War Spatial Autoregressive (SAR)
FML_SAR <-"draw ~ draw.lag + draw.lag.SLfariss + draw.lag.polyarchy + SLreg_civilwar+
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"

model_sar <- milm(fml=FML_SAR, midata=newdata)

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model3)


#Model 4 - Regional Civil War Spatial Lag of X (SLX)

FML_SLX <-"draw ~ draw.lag + draw.lag.SLpolyarchy + draw.lag.polyarchy + SLreg_civilwar+
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_slx <- milm(fml=FML_SLX, midata=newdata)

model4<- data.frame(model_slx$terms, model_slx$beta, model_slx$SE, model_slx$beta/model_slx$SE)
colnames(model4) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model4)


#Model 5 - Regional Civil War Spatial Durbin Model (SDM)
FML_SDM <-"draw ~ draw.lag + draw.lag.SLfariss + draw.lag.SLpolyarchy + draw.lag.polyarchy + SLreg_civilwar+
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"


model_sdm <- milm(fml=FML_SDM, midata=newdata)

model5<- data.frame(model_sdm$terms, model_sdm$beta, model_sdm$SE, model_sdm$beta/model_sdm$SE)
colnames(model5) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model5)

### Table A16 & A17 - Latent measure of conflict ###
### Load data

dat <- read.dta("master_protest.dta")

attach(dat)
datfh <- na.omit(data.frame(theta_mean, theta_sd, l_fariss, l_fariss_sd, SLreg_fariss_t1, SLreg_fariss_sd, v2cltort_sd, v2clkill_sd,
                            l_polity, l_polyarchy, l_polyarchy_sd, SLreg_polyarchy_t1, SLreg_polyarchy_sd, l_uds, l_boix, l_cgv, l_gwf_democracy, 
                            SLreg_protest, SLreg_protest_sd, mean5, sd5,
                            SLreg_polity, SLreg_gwf, SLreg_polyarchy_t1, SLreg_boix_t1, SLreg_cgv_t1,
                            SLreg_clkill_sd, SLreg_clkill_t1, SLreg_cltort_sd, SLreg_cltort_t1, l_v2clkill_sd, l_v2clkill, l_cltort_sd, l_cltort, v2cltort, v2clkill,
                            l_uds2, l_uds_sd, SLreg_uds_t1, SLreg_uds_sd, l_gdp_gro, l_lngdppc, l_lnpop, civil_war, international_war, l_resdep,
                            ccode_fe1, ccode_fe2, ccode_fe3, ccode_fe4, ccode_fe5, ccode_fe6, ccode_fe7, ccode_fe8, ccode_fe9, 
                            ccode_fe10, ccode_fe11, ccode_fe12, ccode_fe13, ccode_fe14, ccode_fe15, ccode_fe16, ccode_fe17, 
                            ccode_fe18, ccode_fe19, ccode_fe20, ccode_fe21, ccode_fe22, ccode_fe23, ccode_fe24, ccode_fe25, 
                            ccode_fe26, ccode_fe27, ccode_fe28, ccode_fe29, ccode_fe30, ccode_fe31, ccode_fe32, ccode_fe33, 
                            ccode_fe34, ccode_fe35, ccode_fe36, ccode_fe37, ccode_fe38, ccode_fe39, ccode_fe40, ccode_fe41, 
                            ccode_fe42, ccode_fe43, ccode_fe44, ccode_fe45, ccode_fe46, ccode_fe47, ccode_fe48, ccode_fe49, 
                            ccode_fe50, ccode_fe51, ccode_fe52, ccode_fe53, ccode_fe54, ccode_fe55, ccode_fe56, ccode_fe57,
                            ccode_fe58, ccode_fe59, ccode_fe60, ccode_fe61, ccode_fe62, ccode_fe63, ccode_fe64, ccode_fe65, 
                            ccode_fe66, ccode_fe67, ccode_fe68, ccode_fe69, ccode_fe70, ccode_fe71, ccode_fe72, ccode_fe73, 
                            ccode_fe74, ccode_fe75, ccode_fe76, ccode_fe77, ccode_fe78, ccode_fe79, ccode_fe80, ccode_fe81, 
                            ccode_fe82, ccode_fe83, ccode_fe84, ccode_fe85, ccode_fe86, ccode_fe87, ccode_fe88, ccode_fe89, 
                            ccode_fe90, ccode_fe91, ccode_fe92, ccode_fe93, ccode_fe94, ccode_fe95, ccode_fe96, ccode_fe97, 
                            ccode_fe98, ccode_fe99, ccode_fe100, ccode_fe101, ccode_fe102, ccode_fe103, ccode_fe104, ccode_fe105, 
                            ccode_fe106, ccode_fe107, ccode_fe108, ccode_fe109, ccode_fe110, ccode_fe111, ccode_fe112, ccode_fe113, 
                            ccode_fe114, ccode_fe115, ccode_fe116, ccode_fe117, ccode_fe118, ccode_fe119, ccode_fe120, ccode_fe121, 
                            ccode_fe122, ccode_fe123, ccode_fe124, ccode_fe125, ccode_fe126, ccode_fe127, ccode_fe128, ccode_fe129, 
                            ccode_fe130, ccode_fe131, ccode_fe132, ccode_fe133, ccode_fe134, ccode_fe135, ccode_fe136, ccode_fe137, 
                            y_fe1, y_fe2, y_fe3, y_fe4, y_fe5, y_fe6, 
                            y_fe7, y_fe8, y_fe9, y_fe10, y_fe11, y_fe12, y_fe13, y_fe14, y_fe15, y_fe16, y_fe17, y_fe18, y_fe19, y_fe20, 
                            y_fe21, y_fe22, y_fe23, y_fe24, y_fe25, y_fe26, y_fe27, y_fe28, y_fe29, y_fe30, y_fe31, y_fe32, y_fe33, y_fe34, 
                            y_fe35, y_fe36, y_fe37, y_fe38, y_fe39, y_fe40, y_fe41, y_fe42, y_fe43, y_fe44, y_fe45, y_fe46, y_fe47, y_fe48, 
                            y_fe49, y_fe50, y_fe51, y_fe52, y_fe53))

dim(datfh)


### Setup data to run models where we account for uncertainty
set.seed(32261991)
newdata <- list()
for(i in 1:1000){
  newdata[[i]] <- datfh
  newdata[[i]]$draw <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$theta_mean, sd=datfh$theta_sd)
  newdata[[i]]$draw.lag <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_fariss, sd=datfh$l_fariss_sd)
  newdata[[i]]$draw.lag.SLfariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_fariss_t1, sd=datfh$SLreg_fariss_sd)
  newdata[[i]]$draw.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2clkill, sd=datfh$v2clkill_sd)
  newdata[[i]]$draw.lag.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_v2clkill, sd=datfh$l_v2clkill_sd)
  newdata[[i]]$draw.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2cltort, sd=datfh$v2cltort_sd)
  newdata[[i]]$draw.lag.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_cltort, sd=datfh$l_cltort_sd)
  newdata[[i]]$draw.lag.SLkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_clkill_t1, sd=datfh$SLreg_clkill_sd)
  newdata[[i]]$draw.lag.SLtort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_cltort_t1, sd=datfh$SLreg_cltort_sd)
  newdata[[i]]$draw.lag.SLpolyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_polyarchy_t1, sd=datfh$SLreg_polyarchy_sd)
  newdata[[i]]$draw.lag.polyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_polyarchy, sd=datfh$l_polyarchy_sd)
  newdata[[i]]$draw.lag.protest <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_protest, sd=datfh$SLreg_protest_sd)
  newdata[[i]]$draw.lag.protest.mean <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$mean5, sd=datfh$sd5)
}


##Model 1 UDS - Non-spatial (NS)
FML_NS <-"draw ~ draw.lag + draw.lag.polyarchy + draw.lag.protest + draw.lag.protest.mean +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep"

model_ns <- milm(fml=FML_NS, midata=newdata)
df <- 5870

# create object for printing during the loop and then print the results to screen
model1<- data.frame(model_ns$terms, model_ns$beta, model_ns$SE, model_ns$beta/model_ns$SE)
colnames(model1) <- c("Variable","Est. Coef.","SE","t-statistic")

print(paste("Model Output"))
print(model1)


#Model 3 - Spatial Autoregressive (SAR)
FML_SAR <-"draw ~ draw.lag + draw.lag.SLfariss + draw.lag.protest + draw.lag.polyarchy + draw.lag.protest.mean +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + 
ccode_fe1+ ccode_fe2+ ccode_fe3+ ccode_fe4+ ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
  ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
  ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
  ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
  ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe40+ ccode_fe41+ 
  ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
  ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
  ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
  ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
  ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
  ccode_fe82+ ccode_fe83+ ccode_fe84+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
  ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
  ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
  ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
  ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
  ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
  ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
  y_fe1+ y_fe3+ y_fe4+ y_fe5+ y_fe6+ 
  y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
  y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
  y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
  y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53"

model_sar <- milm(fml=FML_SAR, midata=newdata)

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model3)
write.table(model3, file = "m3.csv", sep = ",", quote = FALSE, row.names = F)


#Model 4 - Spatial Lag of X (SLX)

FML_SLX <-"draw ~ draw.lag + draw.lag.SLpolyarchy + draw.lag.polyarchy + draw.lag.protest + draw.lag.protest.mean +
l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + 
ccode_fe1+ ccode_fe2+ ccode_fe3+ ccode_fe4+ ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
  ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe40+ ccode_fe41+ 
ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
ccode_fe82+ ccode_fe83+ ccode_fe84+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
y_fe1+ y_fe3+ y_fe4+ y_fe5+ y_fe6+ 
y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53"


model_slx <- milm(fml=FML_SLX, midata=newdata)

model4<- data.frame(model_slx$terms, model_slx$beta, model_slx$SE, model_slx$beta/model_slx$SE)
colnames(model4) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model4)


#Model 5 - Spatial Durbin Model (SDM)
FML_SDM <-"draw ~ draw.lag + draw.lag.SLfariss + draw.lag.SLpolyarchy + draw.lag.polyarchy + draw.lag.protest+
draw.lag.protest.mean +l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + 
ccode_fe1+ ccode_fe2+ ccode_fe3+ ccode_fe4+ ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
  ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
  ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
  ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
  ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe40+ ccode_fe41+ 
  ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
  ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
  ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
  ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
  ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
  ccode_fe82+ ccode_fe83+ ccode_fe84+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
  ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
  ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
  ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
  ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
  ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
  ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
  y_fe1+ y_fe3+ y_fe4+ y_fe5+ y_fe6+ 
  y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
  y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
  y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
  y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53"


model_sdm <- milm(fml=FML_SDM, midata=newdata)

model5<- data.frame(model_sdm$terms, model_sdm$beta, model_sdm$SE, model_sdm$beta/model_sdm$SE)
colnames(model5) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model5)


### Table A19 - Global connectivities m-STAR ###
### Load data

dat <- read.dta("MainData.dta")

attach(dat)
datfh <- na.omit(data.frame(theta_mean, theta_sd, l_fariss, l_fariss_sd, SLreg_fariss_t1, SLreg_fariss_sd, v2cltort_sd, v2clkill_sd,
                            l_polity, l_polyarchy, l_polyarchy_sd, SLreg_polyarchy_t1, SLreg_polyarchy_sd, l_uds, l_boix, l_cgv, l_gwf_democracy, 
                            SLreg_polity, SLreg_gwf, SLreg_polyarchy_t1, SLreg_boix_t1, SLreg_cgv_t1,
                            SLreg_gdp_gro, SLreg_gdppcln, SLreg_lnpop, SLreg_civilwar, SLreg_intwar,
                            SL_def_fariss, SL_def_fariss_sd, SL_dem_fariss, SL_dem_fariss_sd,
                            SLreg_clkill_sd, SLreg_clkill_t1, SLreg_cltort_sd, SLreg_cltort_t1, l_v2clkill_sd, l_v2clkill, l_cltort_sd, l_cltort, v2cltort, v2clkill,
                            l_uds2, l_uds_sd, SLreg_uds_t1, SLreg_uds_sd, l_gdp_gro, l_lngdppc, l_lnpop, civil_war, international_war, l_resdep,
                            ccode_fe1, ccode_fe2, ccode_fe3, ccode_fe4, ccode_fe5, ccode_fe6, ccode_fe7, ccode_fe8, ccode_fe9, 
                            ccode_fe10, ccode_fe11, ccode_fe12, ccode_fe13, ccode_fe14, ccode_fe15, ccode_fe16, ccode_fe17, 
                            ccode_fe18, ccode_fe19, ccode_fe20, ccode_fe21, ccode_fe22, ccode_fe23, ccode_fe24, ccode_fe25, 
                            ccode_fe26, ccode_fe27, ccode_fe28, ccode_fe29, ccode_fe30, ccode_fe31, ccode_fe32, ccode_fe33, 
                            ccode_fe34, ccode_fe35, ccode_fe36, ccode_fe37, ccode_fe38, ccode_fe39, ccode_fe40, ccode_fe41, 
                            ccode_fe42, ccode_fe43, ccode_fe44, ccode_fe45, ccode_fe46, ccode_fe47, ccode_fe48, ccode_fe49, 
                            ccode_fe50, ccode_fe51, ccode_fe52, ccode_fe53, ccode_fe54, ccode_fe55, ccode_fe56, ccode_fe57,
                            ccode_fe58, ccode_fe59, ccode_fe60, ccode_fe61, ccode_fe62, ccode_fe63, ccode_fe64, ccode_fe65, 
                            ccode_fe66, ccode_fe67, ccode_fe68, ccode_fe69, ccode_fe70, ccode_fe71, ccode_fe72, ccode_fe73, 
                            ccode_fe74, ccode_fe75, ccode_fe76, ccode_fe77, ccode_fe78, ccode_fe79, ccode_fe80, ccode_fe81, 
                            ccode_fe82, ccode_fe83, ccode_fe84, ccode_fe85, ccode_fe86, ccode_fe87, ccode_fe88, ccode_fe89, 
                            ccode_fe90, ccode_fe91, ccode_fe92, ccode_fe93, ccode_fe94, ccode_fe95, ccode_fe96, ccode_fe97, 
                            ccode_fe98, ccode_fe99, ccode_fe100, ccode_fe101, ccode_fe102, ccode_fe103, ccode_fe104, ccode_fe105, 
                            ccode_fe106, ccode_fe107, ccode_fe108, ccode_fe109, ccode_fe110, ccode_fe111, ccode_fe112, ccode_fe113, 
                            ccode_fe114, ccode_fe115, ccode_fe116, ccode_fe117, ccode_fe118, ccode_fe119, ccode_fe120, ccode_fe121, 
                            ccode_fe122, ccode_fe123, ccode_fe124, ccode_fe125, ccode_fe126, ccode_fe127, ccode_fe128, ccode_fe129, 
                            ccode_fe130, ccode_fe131, ccode_fe132, ccode_fe133, ccode_fe134, ccode_fe135, ccode_fe136, ccode_fe137, 
                            ccode_fe138, y_fe1, y_fe2, y_fe3, y_fe4, y_fe5, y_fe6, 
                            y_fe7, y_fe8, y_fe9, y_fe10, y_fe11, y_fe12, y_fe13, y_fe14, y_fe15, y_fe16, y_fe17, y_fe18, y_fe19, y_fe20, 
                            y_fe21, y_fe22, y_fe23, y_fe24, y_fe25, y_fe26, y_fe27, y_fe28, y_fe29, y_fe30, y_fe31, y_fe32, y_fe33, y_fe34, 
                            y_fe35, y_fe36, y_fe37, y_fe38, y_fe39, y_fe40, y_fe41, y_fe42, y_fe43, y_fe44, y_fe45, y_fe46, y_fe47, y_fe48, 
                            y_fe49, y_fe50, y_fe51, y_fe52, y_fe53, y_fe54, y_fe55, y_fe56, y_fe57, y_fe58, y_fe59, y_fe60, y_fe61))

dim(datfh)


### Setup data to run models where we account for uncertainty
set.seed(32261991)
newdata <- list()
for(i in 1:1000){
  newdata[[i]] <- datfh
  newdata[[i]]$draw <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$theta_mean, sd=datfh$theta_sd)
  newdata[[i]]$draw.lag <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_fariss, sd=datfh$l_fariss_sd)
  newdata[[i]]$draw.lag.SLfariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_fariss_t1, sd=datfh$SLreg_fariss_sd)
  newdata[[i]]$draw.lag.SLuds <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_uds_t1, sd=datfh$SLreg_uds_sd)
  newdata[[i]]$draw.lag.uds <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_uds2, sd=datfh$l_uds_sd)
  newdata[[i]]$draw.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2clkill, sd=datfh$v2clkill_sd)
  newdata[[i]]$draw.lag.clkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_v2clkill, sd=datfh$l_v2clkill_sd)
  newdata[[i]]$draw.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2cltort, sd=datfh$v2cltort_sd)
  newdata[[i]]$draw.lag.cltort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_cltort, sd=datfh$l_cltort_sd)
  newdata[[i]]$draw.lag.SLkill <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_clkill_t1, sd=datfh$SLreg_clkill_sd)
  newdata[[i]]$draw.lag.SLtort <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_cltort_t1, sd=datfh$SLreg_cltort_sd)
  newdata[[i]]$draw.lag.SLpolyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SLreg_polyarchy_t1, sd=datfh$SLreg_polyarchy_sd)
  newdata[[i]]$draw.lag.polyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_polyarchy, sd=datfh$l_polyarchy_sd)
  newdata[[i]]$draw.lag.SLdem_fariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SL_dem_fariss, sd=datfh$SL_dem_fariss_sd)
  newdata[[i]]$draw.lag.SLdef_fariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$SL_def_fariss, sd=datfh$SL_def_fariss_sd)
}

#Model 3 - Spatial Autoregressive (SAR)
FML_SAR <-"draw ~ draw.lag + draw.lag.SLfariss + draw.lag.SLdem_fariss + draw.lag.SLdef_fariss + 
                      draw.lag.polyarchy + 
                      l_gdp_gro + l_lngdppc + l_lnpop + civil_war + international_war + l_resdep + ccode_fe1+ 
                      ccode_fe2+ ccode_fe3+ ccode_fe4 + ccode_fe5+ ccode_fe6+ ccode_fe7+ ccode_fe8+ ccode_fe9+ 
                      ccode_fe10+ ccode_fe11+ ccode_fe12+ ccode_fe13+ ccode_fe14+ ccode_fe15+ ccode_fe16+ ccode_fe17+ 
                      ccode_fe18+ ccode_fe19+ ccode_fe20+ ccode_fe21+ ccode_fe22+ ccode_fe23+ ccode_fe24+ ccode_fe25+ 
                      ccode_fe26+ ccode_fe27+ ccode_fe28+ ccode_fe29+ ccode_fe30+ ccode_fe31+ ccode_fe32+ ccode_fe33+ 
                      ccode_fe34+ ccode_fe35+ ccode_fe36+ ccode_fe37+ ccode_fe38+ ccode_fe39+ ccode_fe40+ ccode_fe41+ 
                      ccode_fe42+ ccode_fe43+ ccode_fe44+ ccode_fe45+ ccode_fe46+ ccode_fe47+ ccode_fe48+ ccode_fe49+ 
                      ccode_fe50+ ccode_fe51+ ccode_fe52+ ccode_fe53+ ccode_fe54+ ccode_fe55+ ccode_fe56+ ccode_fe57+
                      ccode_fe58+ ccode_fe59+ ccode_fe60+ ccode_fe61+ ccode_fe62+ ccode_fe63+ ccode_fe64+ ccode_fe65+ 
                      ccode_fe66+ ccode_fe67+ ccode_fe68+ ccode_fe69+ ccode_fe70+ ccode_fe71+ ccode_fe72+ ccode_fe73+ 
                      ccode_fe74+ ccode_fe75+ ccode_fe76+ ccode_fe77+ ccode_fe78+ ccode_fe79+ ccode_fe80+ ccode_fe81+ 
                      ccode_fe82+ ccode_fe83+ ccode_fe85+ ccode_fe86+ ccode_fe87+ ccode_fe88+ ccode_fe89+ 
                      ccode_fe90+ ccode_fe91+ ccode_fe92+ ccode_fe93+ ccode_fe94+ ccode_fe95+ ccode_fe96+ ccode_fe97+ 
                      ccode_fe98+ ccode_fe99+ ccode_fe100+ ccode_fe101+ ccode_fe102+ ccode_fe103+ ccode_fe104+ ccode_fe105+ 
                      ccode_fe106+ ccode_fe107+ ccode_fe108+ ccode_fe109+ ccode_fe110+ ccode_fe111+ ccode_fe112+ ccode_fe113+ 
                      ccode_fe114+ ccode_fe115+ ccode_fe116+ ccode_fe117+ ccode_fe118+ ccode_fe119+ ccode_fe120+ ccode_fe121+ 
                      ccode_fe122+ ccode_fe123+ ccode_fe124+ ccode_fe125+ ccode_fe126+ ccode_fe127+ ccode_fe128+ ccode_fe129+ 
                      ccode_fe130+ ccode_fe131+ ccode_fe132+ ccode_fe133+ ccode_fe134+ ccode_fe135+ ccode_fe136+ ccode_fe137+ 
                      ccode_fe138+ y_fe1+ y_fe2+ y_fe3+ y_fe5+ y_fe6+ 
                      y_fe7+ y_fe8+ y_fe9+ y_fe10+ y_fe11+ y_fe12+ y_fe13+ y_fe14+ y_fe15+ y_fe16+ y_fe17+ y_fe18+ y_fe19+ y_fe20+ 
                      y_fe21+ y_fe22+ y_fe23+ y_fe24+ y_fe25+ y_fe26+ y_fe27+ y_fe28+ y_fe29+ y_fe30+ y_fe31+ y_fe32+ y_fe33+ y_fe34+ 
                      y_fe35+ y_fe36+ y_fe37+ y_fe38+ y_fe39+ y_fe40+ y_fe41+ y_fe42+ y_fe43+ y_fe44+ y_fe45+ y_fe46+ y_fe47+ y_fe48+ 
                      y_fe49+ y_fe50+ y_fe51+ y_fe52+ y_fe53+ y_fe54+ y_fe55+ y_fe56+ y_fe57+ y_fe58+ y_fe59+ y_fe60+ y_fe61"

model_sar <- milm(fml=FML_SAR, midata=newdata)

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model3)



### Table A20 - No resources dependence control ###
### Load data
dat <- read.dta("MainData_noresdep.dta")

attach(dat)
datfh <- na.omit(data.frame(v2x_polyarchy, v2x_polyarchy_sd, v2x_clphy, v2x_clphy_sd,
                            l_v2x_polyarchy, l_v2x_polyarchy_sd, l_v2x_clphy, l_v2x_clphy_sd, l_e_migdpgro, 
                            l_e_migdppcln, l_wb_pop1000s, l_wbpopln, l_e_polity2, 
                            l_e_boix_regime, minor, civilwar, int_war, l_minor, l_civilwar, l_int_war, 
                            theta_mean, theta_sd, l_theta_mean, l_theta_sd, uds_mean, uds_sd, l_uds_mean, 
                            l_uds_sd, gwf_democracy, l_gwf_dem, slreg_polyarchy_t1, slreg_polyarchy_sd_t1, 
                            slreg_fariss_t1, slreg_fariss_sd_t1, slreg_clphy_t1, slreg_clphy_sd_t1, slreg_uds_t1, 
                            slreg_uds_sd_t1, slreg_polity_t1, slreg_bmr_t1, slreg_gwf_t1, slmind_polyarchy_t1, 
                            slmind_polyarchy_sd_t1, slmind_fariss_t1, slmind_fariss_sd_t1, slmind_clphy_t1, 
                            slmind_clphy_sd_t1, slmind_uds_t1, slmind_uds_sd_t1, slmind_polity_t1, slmind_gwf_t1,
                            ccode_fe1 ,ccode_fe2 ,ccode_fe3 ,ccode_fe4 ,ccode_fe5 ,ccode_fe6 ,ccode_fe7, 
                            ccode_fe8 ,ccode_fe9 ,ccode_fe10 ,ccode_fe11 ,ccode_fe12 ,ccode_fe13 ,ccode_fe14, 
                            ccode_fe15 ,ccode_fe16 ,ccode_fe17 ,ccode_fe18 ,ccode_fe19 ,ccode_fe20 ,ccode_fe21, 
                            ccode_fe22 ,ccode_fe23 ,ccode_fe24 ,ccode_fe25 ,ccode_fe26 ,ccode_fe27 ,ccode_fe28, 
                            ccode_fe29 ,ccode_fe30 ,ccode_fe31 ,ccode_fe32 ,ccode_fe33 ,ccode_fe34 ,ccode_fe35, 
                            ccode_fe36 ,ccode_fe37 ,ccode_fe38 ,ccode_fe39 ,ccode_fe40 ,ccode_fe41 ,ccode_fe42, 
                            ccode_fe43 ,ccode_fe44 ,ccode_fe45 ,ccode_fe46 ,ccode_fe47 ,ccode_fe48 ,ccode_fe49, 
                            ccode_fe50 ,ccode_fe51 ,ccode_fe52 ,ccode_fe53 ,ccode_fe54 ,ccode_fe55 ,ccode_fe56, 
                            ccode_fe57 ,ccode_fe58 ,ccode_fe59 ,ccode_fe60 ,ccode_fe61 ,ccode_fe62 ,ccode_fe63, 
                            ccode_fe64 ,ccode_fe65 ,ccode_fe66 ,ccode_fe67 ,ccode_fe68 ,ccode_fe69 ,ccode_fe70, 
                            ccode_fe71 ,ccode_fe72 ,ccode_fe73 ,ccode_fe74 ,ccode_fe75 ,ccode_fe76 ,ccode_fe77, 
                            ccode_fe78 ,ccode_fe79 ,ccode_fe80 ,ccode_fe81 ,ccode_fe82 ,ccode_fe83 ,ccode_fe84, 
                            ccode_fe85 ,ccode_fe86 ,ccode_fe87 ,ccode_fe88 ,ccode_fe89 ,ccode_fe90 ,ccode_fe91, 
                            ccode_fe92 ,ccode_fe93 ,ccode_fe94 ,ccode_fe95 ,ccode_fe96 ,ccode_fe97 ,ccode_fe98, 
                            ccode_fe99 ,ccode_fe100 ,ccode_fe101 ,ccode_fe102 ,ccode_fe103 ,ccode_fe104 ,ccode_fe105, 
                            ccode_fe106 ,ccode_fe107 ,ccode_fe108 ,ccode_fe109 ,ccode_fe110 ,ccode_fe111, 
                            ccode_fe112 ,ccode_fe113 ,ccode_fe114 ,ccode_fe115 ,ccode_fe116 ,ccode_fe117, 
                            ccode_fe118 ,ccode_fe119 ,ccode_fe120 ,ccode_fe121 ,ccode_fe122 ,ccode_fe123, 
                            ccode_fe124 ,ccode_fe125 ,ccode_fe126 ,ccode_fe127 ,ccode_fe128 ,ccode_fe129, 
                            ccode_fe130 ,ccode_fe131 ,ccode_fe132 ,ccode_fe133 ,ccode_fe134 ,ccode_fe135, 
                            ccode_fe136 ,ccode_fe137 ,ccode_fe138 ,ccode_fe139 ,ccode_fe140 ,ccode_fe141, 
                            ccode_fe142 ,ccode_fe143 ,ccode_fe144 ,ccode_fe145 ,ccode_fe146 ,ccode_fe147, 
                            ccode_fe148 ,ccode_fe149 ,ccode_fe150 ,ccode_fe151 ,ccode_fe152 ,ccode_fe153, 
                            ccode_fe154 ,ccode_fe155 ,ccode_fe156 ,ccode_fe157 ,ccode_fe158 ,ccode_fe159, 
                            ccode_fe160 ,ccode_fe161 ,ccode_fe162 ,ccode_fe163,
                            ccode_fe164 ,ccode_fe165 ,ccode_fe166 ,ccode_fe167, ccode_fe168,
                            y_fe1 ,y_fe2 ,y_fe3 ,y_fe4 ,y_fe5 ,y_fe6 ,y_fe7 ,y_fe8 ,y_fe9 ,y_fe10 ,y_fe11 ,
                            y_fe12 ,y_fe13 ,y_fe14 ,y_fe15 ,y_fe16 ,y_fe17 ,y_fe18 ,y_fe19 ,y_fe20 ,y_fe21 ,
                            y_fe22 ,y_fe23 ,y_fe24 ,y_fe25 ,y_fe26 ,y_fe27 ,y_fe28 ,y_fe29 ,y_fe30 ,y_fe31 ,
                            y_fe32 ,y_fe33 ,y_fe34 ,y_fe35 ,y_fe36 ,y_fe37 ,y_fe38 ,y_fe39 ,y_fe40 ,y_fe41 ,
                            y_fe42 ,y_fe43 ,y_fe44 ,y_fe45 ,y_fe46 ,y_fe47 ,y_fe48 ,y_fe49 ,y_fe50 ,y_fe51 ,
                            y_fe52 ,y_fe53 ,y_fe54 ,y_fe55 ,y_fe56 ,y_fe57 ,y_fe58 ,y_fe59 ,y_fe60 ,y_fe61))

dim(datfh)


### Setup data to run models where we account for uncertainty
set.seed(32261991)
newdata <- list()
for(i in 1:1000){
  newdata[[i]] <- datfh
  newdata[[i]]$draw <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$theta_mean, sd=datfh$theta_sd)
  newdata[[i]]$draw.lag <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_theta_mean, sd=datfh$l_theta_sd)
  newdata[[i]]$draw.lag.SLfariss <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$slreg_fariss_t1, sd=datfh$slreg_fariss_sd_t1)
  newdata[[i]]$draw.clphy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2x_clphy, sd=datfh$v2x_clphy_sd)
  newdata[[i]]$draw.lag.clphy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_v2x_clphy, sd=datfh$l_v2x_clphy_sd)
  newdata[[i]]$draw.lag.SLclphy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$slreg_clphy_t1, sd=datfh$slreg_clphy_sd_t1)
  newdata[[i]]$draw.polyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$v2x_polyarchy, sd=datfh$v2x_polyarchy_sd)
  newdata[[i]]$draw.lag.SLpolyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$slreg_polyarchy_t1, sd=datfh$slreg_polyarchy_sd_t1)
  newdata[[i]]$draw.lag.polyarchy <- rnorm(cbind(rep(1,nrow(datfh))), mean=datfh$l_v2x_polyarchy, sd=datfh$l_v2x_polyarchy_sd)
}



#################### Polyarchy ##########################

#Model 1 - NS 
FML_NS <- "draw ~ draw.lag + draw.lag.polyarchy + 
                  l_e_migdpgro + l_e_migdppcln + l_wbpopln +
                  civilwar + int_war"

model_ns <- milm(fml=FML_NS, midata=newdata)
df <- 7137

model1<- data.frame(model_ns$terms, model_ns$beta, model_ns$SE, model_ns$beta/model_ns$SE)
colnames(model1) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model1)


#Model 2 - NS + FE
FML_NS_FE <-"draw ~ draw.lag + draw.lag.polyarchy + 
                  l_e_migdpgro + l_e_migdppcln + l_wbpopln +
                  civilwar + int_war+ 
      ccode_fe1 +ccode_fe2 +ccode_fe3 +ccode_fe4 +ccode_fe5 +ccode_fe6 +ccode_fe7+ 
      ccode_fe8 +ccode_fe9 +ccode_fe10 +ccode_fe11 +ccode_fe12 +ccode_fe13 +ccode_fe14 + 
      ccode_fe15 +ccode_fe16 +ccode_fe19 +ccode_fe20 +ccode_fe21 +
      ccode_fe22 +ccode_fe23 +ccode_fe24 +ccode_fe25 +ccode_fe26 +ccode_fe27 +ccode_fe28 +
      ccode_fe29 +ccode_fe30 +ccode_fe31 +ccode_fe32 +ccode_fe33 +ccode_fe34 +ccode_fe35 +
      ccode_fe36 +ccode_fe37 +ccode_fe38 +ccode_fe39 +ccode_fe40 +ccode_fe41 +ccode_fe42 +
      ccode_fe43 +ccode_fe44 +ccode_fe45 +ccode_fe46 +ccode_fe48 +ccode_fe49 +
      ccode_fe50 +ccode_fe51 +ccode_fe52 +ccode_fe53 +ccode_fe54 +ccode_fe55 +ccode_fe56 +
      ccode_fe57 +ccode_fe58 +ccode_fe59 +ccode_fe60 +ccode_fe61 +ccode_fe62 +ccode_fe63 +
      ccode_fe64 +ccode_fe65 +ccode_fe66 +ccode_fe67 +ccode_fe68 +ccode_fe69 +ccode_fe70 +
      ccode_fe71 +ccode_fe72 +ccode_fe73 +ccode_fe74 +ccode_fe75 +ccode_fe76 +ccode_fe77 +
      ccode_fe78 +ccode_fe79 +ccode_fe80 +ccode_fe81 +ccode_fe82 +ccode_fe83 +ccode_fe84 +
      ccode_fe85 +ccode_fe86 +ccode_fe87 +ccode_fe88 +ccode_fe89 +ccode_fe90 +ccode_fe91 +
      ccode_fe92 +ccode_fe93 +ccode_fe94 +ccode_fe96 +ccode_fe97  +
      ccode_fe99 +ccode_fe100 +ccode_fe101 +ccode_fe102 +ccode_fe103 +ccode_fe104 +ccode_fe105 +
      ccode_fe106 +ccode_fe107 +ccode_fe108 +ccode_fe109 +ccode_fe110 +ccode_fe111 +
      ccode_fe112 +ccode_fe113 +ccode_fe114 +ccode_fe115 +ccode_fe116 +ccode_fe117 +
      ccode_fe118 +ccode_fe119 +ccode_fe120 +ccode_fe121 +ccode_fe122 +ccode_fe123 +
      ccode_fe124 +ccode_fe125 +ccode_fe126 +ccode_fe127 +ccode_fe128 +ccode_fe129 +
      ccode_fe130 +ccode_fe131 +ccode_fe132 +ccode_fe133 +ccode_fe134 +ccode_fe135 +
      ccode_fe136 +ccode_fe137 +ccode_fe138 +ccode_fe139 +ccode_fe140 +ccode_fe141 +
      ccode_fe142 +ccode_fe143 +ccode_fe145 +ccode_fe146 +ccode_fe147 +
      ccode_fe148 +ccode_fe149 +ccode_fe150 +ccode_fe151 +ccode_fe152 +ccode_fe153 + 
      ccode_fe154 +ccode_fe155 +ccode_fe156 +ccode_fe157 +ccode_fe159 +
      ccode_fe161 + y_fe1 +y_fe3 +y_fe4 +y_fe5 +y_fe6 +y_fe7 +y_fe8 +y_fe9 +y_fe10 +y_fe11 + 
    y_fe12 +y_fe13 +y_fe14 +y_fe15 +y_fe16 +y_fe17 +y_fe18 +y_fe19 +y_fe20 +y_fe21 +
    y_fe22 +y_fe23 +y_fe24 +y_fe25 +y_fe26 +y_fe27 +y_fe28 +y_fe29 +y_fe30 +y_fe31 +
    y_fe32 +y_fe33 +y_fe34 +y_fe35 +y_fe36 +y_fe37 +y_fe38 +y_fe39 +y_fe40 +y_fe41 +
    y_fe42 +y_fe43 +y_fe44 +y_fe45 +y_fe46 +y_fe47 +y_fe48 +y_fe49 +y_fe50 +y_fe51 +
    y_fe52 +y_fe53 +y_fe54 +y_fe55 +y_fe56 +y_fe57 +y_fe58 +y_fe59 +y_fe60 +y_fe61"

model_ns_fe <- milm(fml=FML_NS_FE, midata=newdata)
df <- 7137

model2<- data.frame(model_ns_fe$terms, model_ns_fe$beta, model_ns_fe$SE, model_ns_fe$beta/model_ns_fe$SE)
colnames(model2) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model2)

#Model 3 - Spatial Autoregressive (SAR)
FML_SAR <-"draw ~ draw.lag + draw.lag.polyarchy + draw.lag.SLfariss +
                  l_e_migdpgro + l_e_migdppcln + l_wbpopln +
                  civilwar + int_war+ 
      ccode_fe1 +ccode_fe2 +ccode_fe3 +ccode_fe4 +ccode_fe5 +ccode_fe6 +ccode_fe7+ 
      ccode_fe8 +ccode_fe9 +ccode_fe10 +ccode_fe11 +ccode_fe12 +ccode_fe13 +ccode_fe14 + 
      ccode_fe15 +ccode_fe16 +ccode_fe19 +ccode_fe20 +ccode_fe21 +
      ccode_fe22 +ccode_fe23 +ccode_fe24 +ccode_fe25 +ccode_fe26 +ccode_fe27 +ccode_fe28 +
      ccode_fe29 +ccode_fe30 +ccode_fe31 +ccode_fe32 +ccode_fe33 +ccode_fe34 +ccode_fe35 +
      ccode_fe36 +ccode_fe37 +ccode_fe38 +ccode_fe39 +ccode_fe40 +ccode_fe41 +ccode_fe42 +
      ccode_fe43 +ccode_fe44 +ccode_fe45 +ccode_fe46 +ccode_fe48 +ccode_fe49 +
      ccode_fe50 +ccode_fe51 +ccode_fe52 +ccode_fe53 +ccode_fe54 +ccode_fe55 +ccode_fe56 +
      ccode_fe57 +ccode_fe58 +ccode_fe59 +ccode_fe60 +ccode_fe61 +ccode_fe62 +ccode_fe63 +
      ccode_fe64 +ccode_fe65 +ccode_fe66 +ccode_fe67 +ccode_fe68 +ccode_fe69 +ccode_fe70 +
      ccode_fe71 +ccode_fe72 +ccode_fe73 +ccode_fe74 +ccode_fe75 +ccode_fe76 +ccode_fe77 +
      ccode_fe78 +ccode_fe79 +ccode_fe80 +ccode_fe81 +ccode_fe82 +ccode_fe83 +ccode_fe84 +
      ccode_fe85 +ccode_fe86 +ccode_fe87 +ccode_fe88 +ccode_fe89 +ccode_fe90 +ccode_fe91 +
      ccode_fe92 +ccode_fe93 +ccode_fe94 +ccode_fe96 +ccode_fe97  +
      ccode_fe99 +ccode_fe100 +ccode_fe101 +ccode_fe102 +ccode_fe103 +ccode_fe104 +ccode_fe105 +
      ccode_fe106 +ccode_fe107 +ccode_fe108 +ccode_fe109 +ccode_fe110 +ccode_fe111 +
      ccode_fe112 +ccode_fe113 +ccode_fe114 +ccode_fe115 +ccode_fe116 +ccode_fe117 +
      ccode_fe118 +ccode_fe119 +ccode_fe120 +ccode_fe121 +ccode_fe122 +ccode_fe123 +
      ccode_fe124 +ccode_fe125 +ccode_fe126 +ccode_fe127 +ccode_fe128 +ccode_fe129 +
      ccode_fe130 +ccode_fe131 +ccode_fe132 +ccode_fe133 +ccode_fe134 +ccode_fe135 +
      ccode_fe136 +ccode_fe137 +ccode_fe138 +ccode_fe139 +ccode_fe140 +ccode_fe141 +
      ccode_fe142 +ccode_fe143 +ccode_fe145 +ccode_fe146 +ccode_fe147 +
      ccode_fe148 +ccode_fe149 +ccode_fe150 +ccode_fe151 +ccode_fe152 +ccode_fe153 + 
      ccode_fe154 +ccode_fe155 +ccode_fe156 +ccode_fe157 +ccode_fe159 +
      ccode_fe161 + y_fe1 +y_fe3 +y_fe4 +y_fe5 +y_fe6 +y_fe7 +y_fe8 +y_fe9 +y_fe10 +y_fe11 + 
    y_fe12 +y_fe13 +y_fe14 +y_fe15 +y_fe16 +y_fe17 +y_fe18 +y_fe19 +y_fe20 +y_fe21 +
    y_fe22 +y_fe23 +y_fe24 +y_fe25 +y_fe26 +y_fe27 +y_fe28 +y_fe29 +y_fe30 +y_fe31 +
    y_fe32 +y_fe33 +y_fe34 +y_fe35 +y_fe36 +y_fe37 +y_fe38 +y_fe39 +y_fe40 +y_fe41 +
    y_fe42 +y_fe43 +y_fe44 +y_fe45 +y_fe46 +y_fe47 +y_fe48 +y_fe49 +y_fe50 +y_fe51 +
    y_fe52 +y_fe53 +y_fe54 +y_fe55 +y_fe56 +y_fe57 +y_fe58 +y_fe59 +y_fe60 +y_fe61"

model_sar <- milm(fml=FML_SAR, midata=newdata)
df <- 7137

model3<- data.frame(model_sar$terms, model_sar$beta, model_sar$SE, model_sar$beta/model_sar$SE)
colnames(model3) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model3)

#Model 4 - Spatial Lag of X (SLX)
FML_SLX <-"draw ~ draw.lag + draw.lag.polyarchy + draw.lag.SLpolyarchy +
                  l_e_migdpgro + l_e_migdppcln + l_wbpopln +
                  civilwar + int_war+ 
      ccode_fe1 +ccode_fe2 +ccode_fe3 +ccode_fe4 +ccode_fe5 +ccode_fe6 +ccode_fe7+ 
      ccode_fe8 +ccode_fe9 +ccode_fe10 +ccode_fe11 +ccode_fe12 +ccode_fe13 +ccode_fe14 + 
      ccode_fe15 +ccode_fe16 +ccode_fe19 +ccode_fe20 +ccode_fe21 +
      ccode_fe22 +ccode_fe23 +ccode_fe24 +ccode_fe25 +ccode_fe26 +ccode_fe27 +ccode_fe28 +
      ccode_fe29 +ccode_fe30 +ccode_fe31 +ccode_fe32 +ccode_fe33 +ccode_fe34 +ccode_fe35 +
      ccode_fe36 +ccode_fe37 +ccode_fe38 +ccode_fe39 +ccode_fe40 +ccode_fe41 +ccode_fe42 +
      ccode_fe43 +ccode_fe44 +ccode_fe45 +ccode_fe46 +ccode_fe48 +ccode_fe49 +
      ccode_fe50 +ccode_fe51 +ccode_fe52 +ccode_fe53 +ccode_fe54 +ccode_fe55 +ccode_fe56 +
      ccode_fe57 +ccode_fe58 +ccode_fe59 +ccode_fe60 +ccode_fe61 +ccode_fe62 +ccode_fe63 +
      ccode_fe64 +ccode_fe65 +ccode_fe66 +ccode_fe67 +ccode_fe68 +ccode_fe69 +ccode_fe70 +
      ccode_fe71 +ccode_fe72 +ccode_fe73 +ccode_fe74 +ccode_fe75 +ccode_fe76 +ccode_fe77 +
      ccode_fe78 +ccode_fe79 +ccode_fe80 +ccode_fe81 +ccode_fe82 +ccode_fe83 +ccode_fe84 +
      ccode_fe85 +ccode_fe86 +ccode_fe87 +ccode_fe88 +ccode_fe89 +ccode_fe90 +ccode_fe91 +
      ccode_fe92 +ccode_fe93 +ccode_fe94 +ccode_fe96 +ccode_fe97  +
      ccode_fe99 +ccode_fe100 +ccode_fe101 +ccode_fe102 +ccode_fe103 +ccode_fe104 +ccode_fe105 +
      ccode_fe106 +ccode_fe107 +ccode_fe108 +ccode_fe109 +ccode_fe110 +ccode_fe111 +
      ccode_fe112 +ccode_fe113 +ccode_fe114 +ccode_fe115 +ccode_fe116 +ccode_fe117 +
      ccode_fe118 +ccode_fe119 +ccode_fe120 +ccode_fe121 +ccode_fe122 +ccode_fe123 +
      ccode_fe124 +ccode_fe125 +ccode_fe126 +ccode_fe127 +ccode_fe128 +ccode_fe129 +
      ccode_fe130 +ccode_fe131 +ccode_fe132 +ccode_fe133 +ccode_fe134 +ccode_fe135 +
      ccode_fe136 +ccode_fe137 +ccode_fe138 +ccode_fe139 +ccode_fe140 +ccode_fe141 +
      ccode_fe142 +ccode_fe143 +ccode_fe145 +ccode_fe146 +ccode_fe147 +
      ccode_fe148 +ccode_fe149 +ccode_fe150 +ccode_fe151 +ccode_fe152 +ccode_fe153 + 
      ccode_fe154 +ccode_fe155 +ccode_fe156 +ccode_fe157 +ccode_fe159 +
      ccode_fe161 + y_fe1 +y_fe3 +y_fe4 +y_fe5 +y_fe6 +y_fe7 +y_fe8 +y_fe9 +y_fe10 +y_fe11 + 
    y_fe12 +y_fe13 +y_fe14 +y_fe15 +y_fe16 +y_fe17 +y_fe18 +y_fe19 +y_fe20 +y_fe21 +
    y_fe22 +y_fe23 +y_fe24 +y_fe25 +y_fe26 +y_fe27 +y_fe28 +y_fe29 +y_fe30 +y_fe31 +
    y_fe32 +y_fe33 +y_fe34 +y_fe35 +y_fe36 +y_fe37 +y_fe38 +y_fe39 +y_fe40 +y_fe41 +
    y_fe42 +y_fe43 +y_fe44 +y_fe45 +y_fe46 +y_fe47 +y_fe48 +y_fe49 +y_fe50 +y_fe51 +
    y_fe52 +y_fe53 +y_fe54 +y_fe55 +y_fe56 +y_fe57 +y_fe58 +y_fe59 +y_fe60 +y_fe61"

model_slx <- milm(fml=FML_SLX, midata=newdata)
df <- 7137

model4<- data.frame(model_slx$terms, model_slx$beta, model_slx$SE, model_slx$beta/model_slx$SE)
colnames(model4) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model4)


#Model 5 - Spatial Durbin Model (SDM)
FML_SDM <- "draw ~ draw.lag + draw.lag.polyarchy + draw.lag.SLfariss +
                  draw.lag.SLpolyarchy + l_e_migdpgro + l_e_migdppcln + l_wbpopln +
                  civilwar + int_war+ 
      ccode_fe1 +ccode_fe2 +ccode_fe3 +ccode_fe4 +ccode_fe5 +ccode_fe6 +ccode_fe7+ 
      ccode_fe8 +ccode_fe9 +ccode_fe10 +ccode_fe11 +ccode_fe12 +ccode_fe13 +ccode_fe14 + 
      ccode_fe15 +ccode_fe16 +ccode_fe19 +ccode_fe20 +ccode_fe21 +
      ccode_fe22 +ccode_fe23 +ccode_fe24 +ccode_fe25 +ccode_fe26 +ccode_fe27 +ccode_fe28 +
      ccode_fe29 +ccode_fe30 +ccode_fe31 +ccode_fe32 +ccode_fe33 +ccode_fe34 +ccode_fe35 +
      ccode_fe36 +ccode_fe37 +ccode_fe38 +ccode_fe39 +ccode_fe40 +ccode_fe41 +ccode_fe42 +
      ccode_fe47 +ccode_fe44 +ccode_fe45 +ccode_fe46 +ccode_fe48 +ccode_fe49 +
      ccode_fe50 +ccode_fe51 +ccode_fe52 +ccode_fe53 +ccode_fe54 +ccode_fe55 +ccode_fe56 +
      ccode_fe57 +ccode_fe58 +ccode_fe59 +ccode_fe60 +ccode_fe61 +ccode_fe62 +ccode_fe63 +
      ccode_fe64 +ccode_fe65 +ccode_fe66 +ccode_fe67 +ccode_fe68 +ccode_fe69 +ccode_fe70 +
      ccode_fe71 +ccode_fe72 +ccode_fe73 +ccode_fe74 +ccode_fe75 +ccode_fe76 +ccode_fe77 +
      ccode_fe78 +ccode_fe79 +ccode_fe80 +ccode_fe81 +ccode_fe82 +ccode_fe83 +ccode_fe84 +
      ccode_fe85 +ccode_fe86 +ccode_fe87 +ccode_fe88 +ccode_fe89 +ccode_fe90 +ccode_fe91 +
      ccode_fe92 +ccode_fe93 +ccode_fe94 +ccode_fe96 +ccode_fe97  +
      ccode_fe99 +ccode_fe100 +ccode_fe101 +ccode_fe102 +ccode_fe103 +ccode_fe104 +ccode_fe105 +
      ccode_fe106 +ccode_fe107 +ccode_fe108 +ccode_fe109 +ccode_fe110 +ccode_fe111 +
      ccode_fe112 +ccode_fe113 +ccode_fe114 +ccode_fe115 +ccode_fe116 +ccode_fe117 +
      ccode_fe118 +ccode_fe119 +ccode_fe120 +ccode_fe121 +ccode_fe122 +ccode_fe123 +
      ccode_fe124 +ccode_fe125 +ccode_fe126 +ccode_fe127 +ccode_fe128 +ccode_fe129 +
      ccode_fe130 +ccode_fe131 +ccode_fe132 +ccode_fe133 +ccode_fe134 +ccode_fe135 +
      ccode_fe136 +ccode_fe137 +ccode_fe138 +ccode_fe139 +ccode_fe140 +ccode_fe141 +
      ccode_fe142 +ccode_fe143 +ccode_fe145 +ccode_fe146 +ccode_fe147 +
      ccode_fe148 +ccode_fe149 +ccode_fe150 +ccode_fe151 +ccode_fe152 +ccode_fe153 + 
      ccode_fe154 +ccode_fe155 +ccode_fe156 +ccode_fe157 +ccode_fe159 +
      ccode_fe161 + y_fe1 +y_fe3 +y_fe4 +y_fe5 +y_fe6 +y_fe7 +y_fe8 +y_fe9 +y_fe10 +y_fe11 + 
    y_fe12 +y_fe13 +y_fe14 +y_fe15 +y_fe16 +y_fe17 +y_fe18 +y_fe19 +y_fe20 +y_fe21 +
    y_fe22 +y_fe23 +y_fe24 +y_fe25 +y_fe26 +y_fe27 +y_fe28 +y_fe29 +y_fe30 +y_fe31 +
    y_fe32 +y_fe33 +y_fe34 +y_fe35 +y_fe36 +y_fe37 +y_fe38 +y_fe39 +y_fe40 +y_fe41 +
    y_fe42 +y_fe43 +y_fe44 +y_fe45 +y_fe46 +y_fe47 +y_fe48 +y_fe49 +y_fe50 +y_fe51 +
    y_fe52 +y_fe53 +y_fe54 +y_fe55 +y_fe56 +y_fe57 +y_fe58 +y_fe59 +y_fe60 +y_fe61"

model_sdm <- milm(fml=FML_SDM, midata=newdata)
df<-7137

model5<- data.frame(model_sdm$terms, model_sdm$beta, model_sdm$SE, model_sdm$beta/model_sdm$SE)
colnames(model5) <- c("Variable","Est. Coef.","SE","z-score")
print(paste("Model Output"))
print(model5)
